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Abstract 



We introduce multiple Wilson polynomials, which give a new example of multiple orthogonal 



polynomials (Hermite-Pade polynomials) of type II. These polynomials can be written as 

■ a Jacobi-Piheiro transform, which is a generalization of the Jacobi transform for Wilson 
| polynomials, found by T.H. Koornwinder. Here we need to introduce Jacobi and Jacobi- 

■ Pineiro polynomials with complex parameters. Some explicit formulas are provided for both 
Jacobi-Piheiro and multiple Wilson polynomials, one of them in terms of Kampe de Feriet 
series. Finally we look at some limiting relations and construct a part of a multiple AT-Askey 
table. 



1 Introduction 



Multiple orthogonal polynomials are a generalization of orthogonal polynomials in the sense that 
they satisfy orthogonality conditions with respect to m £ N measures fj,x, ■ ■ ■ ,/U m [22, 3, 14]. In 
this paper m will always represent the number of weights. Multiple orthogonal polynomials arise 
naturally in the theory of simultaneous rational approximation, in particular in Hermite-Pade 
approximation of a system of m (Markov) functions [20, 7, 8]. 

There are two types of multiple orthogonal polynomials. In the present paper we only con- 
sider multiple orthogonal polynomials of type II. Let n = (m, ri2, ■ ■ ■ ,n m ) be a vector of m 
nonnegative integers, which is called a multi-index with length \n\ := n\ + n<i + • • • + n m . Fur- 
thermore let Ti, . . . , T m be the supports of the m measures. A multiple orthogonal polynomial 
Pfl of type II with respect to the multi-index n, is a (nontrivial) polynomial of degree < |n| 
which satisfies the orthogonality conditions 

PH (t)t k d H (t) = 0, k = 0,..., nj -l, j = l,...,m. (1.1) 
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Notice that the measures in (1.1) are not necessarily supposed to be positive. In case of a complex 
orthogonality relation, one usually refers to J)fi clS cl formal multiple orthogonal polynomial. 

Equation (1.1) leads to a system of |n| homogeneous linear relations for the |n| + 1 unknown 
coefficients of p^. A basic requirement in the study of such multiple orthogonal polynomials 
is that there is (up to a scalar multiplicative constant) a unique solution of system (1.1). We 
call n a normal index for ^i, ... , fi m if any solution of (1.1) has exact degree \n\ (which implies 
uniqueness). Let mj^ = J T . t k d[Xj be the kth moment of the measure [Xj. Further set 
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is an nj x (|n|) matrix of moments of the measure fij. Then is the matrix of the linear 
system (1.1) without the last column. It is known and easily verified that the multi-index 
n = (n±, . . . ,n r ) is normal if and only if this matrix has rank \n\. When every multi-index 
is normal we call the system of measures a perfect system. For perfect systems, the multiple 
orthogonal polynomials of type II satisfy a recurrence relation of order m+ 1. The proof is similar 
as in the case of orthogonal polynomials, see for instance [3]. Because of this recurrence relation, 
formal multiple orthogonal polynomials are a useful tool in the spectral theory of non-symmetric 
linear difference operators. 

In the literature one can find some examples of multiple orthogonal polynomials with respect 
to positive measures on the real line which have the same flavor as the classical orthogonal 
polynomials. Two classes of measures have been analyzed in more detail and are known to form 
a perfect system, see for instance the monograph [22] or the survey given in [14]. The first 
class consists of Angelesco systems where the supports of the measures are disjoint intervals. 
In the second class of so-called AT systems, the supports of the m measures coincide, and the 
functions djij{x) / d^i\{x) for j = 1, . . . ,m form an algebraic Chebyshev system [22, Section IV. 4] 
on the convex hull of the support. In the continuous case (where the measures can be written as 
d/j,j(x) = Wj(x)dx, with Wj the weight function of the measure jij) there are multiple Hermite, 
multiple Laguerre I and II, Jacobi-Piheiro, multiple Bessel, Jacobi-Angeleso, Jacobi-Laguerre 
and Laguerre-Hermite polynomials, see [22, 14, 4] and the references therein. Some classical 
discrete examples are multiple Charlier, multiple Kravchuk, multiple Meixner I and II and 
multiple Hahn [5] . For all these examples there exists a first-order raising operator, based on the 
differential operator D or the difference operators A and V, and a Rodrigues formula so that they 
can be called classical. Moreover, there exist differential or difference equations of order m + 1 
(with polynomial coefficients) [4] . The recurrence relations of order m + 1 are known explicitly 
for these examples in the case m < 2. Finally we mention that there also exist some examples 
of multiple orthogonal polynomials associated with modified Bessel functions [12, 13, 26] which 
can be called classical. 

In Subsection 3.1 we recall the definition of one of these examples, namely Jacobi-Piheiro 



polynomials P~*'^\ which are orthogonal with respect to the weights Wj(x) 



'(1 



x) 13 on 



[0,1], aj,P > —1. These polynomials reduce to the classical Jacobi polynomials (shifted to the 
interval [0,1]) for m = 1. We show in Subsection 2.1 that Jacobi polynomials remain formal 
orthogonal polynomials for complex parameters (Xi,(3, the corresponding complex orthogonality 
relation being obtained via an analytic extension of the Beta function in both variables. As we 
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show in Subsection 3.1, also Jacobi-Piheiro polynomials with complex parameters are formal 
multiple orthogonal polynomials of type II. 

In Subsection 3.2 we then introduce the formal multiple Wilson polynomials pft(-;a,b,c,d), 
which give a new example of formal multiple orthogonal polynomials of type II. They are an 
extension of the formal Wilson polynomials p n (-; a, b, c, d) [27, 28] for which we recall the defini- 
tion in Subsection 2.2. We also mention that, with some conditions on the complex parameters 
a,b,c, d, we find the Wilson and Racah polynomials on the top of the Askey scheme which have 
real orthogonality conditions. 

The formal multiple Wilson polynomials satisfy a complex orthogonality conditions with 
respect to m Wilson weights 

Tja + t)T(a - t)T(b 3 + QTfo - t)T(c + t)T(c - t)T(d + t)T(d - t) 
w[t,a,0j,c,a) - r(2i)r(-2i) ' [ j 

j = 1, . . . , m, where we integrate over the imaginary axis deformed so as to separate the increas- 
ing sequences of poles of these weight functions from the decreasing ones. Note that there are 
some additional conditions on the complex parameters in order to have that these Wilson weights 
have only simple poles. We prove in Theorem 3.3 that the weight functions (1.3) form a perfect 
system if bi — bt ^ Z whenever I ^ k. In the same Theorem we show that, for 3ft(c + d) > and 
< |5ft(i)| < 3i(a), the formal multiple Wilson polynomials can be written as the Jacobi-Piheiro 
transform 

p H (t 2 ;a,b,c,d)=k H P^' (3) (u)w {a - 1 ^(u)K(u,t;a,0,c,d)du, (1.4) 
Jo 

where a = (a + b\ — 1, a + b m — 1) and (3 = c + d — 1. Here k^ is a normalizing constant, 
w( a 'P\u) = u a (l — u)^ the Jacobi weight and 



u- h - 1 „ ( c-t,d-t 



1-u) (1.5) 



K(u, t; a, b, c, d) = — — — — , —. , 2*1 . 

F(a — t)F(a + t)T(c + d) \ c + d 

some kernel function, independent of n. In the scalar case (m = 1) this formula reduces to a 
Jacobi transform for the Wilson polynomials which was already found by T.H. Koornwinder [18, 
(3.3)]. We recall this Jacobi transform in Subsection 2.3 and give a short proof. Note that the 
parameters a,b±, . . . , b m , c, d can take complex values. 

The Jacobi-Piheiro transform (1.4) is the key formula of this paper. In Subsection 3.1 
we obtain two new hypergeometric representations for the Jacobi-Piheiro polynomials starting 
from the Rodrigues formula. Applying the Jacobi-Piheiro transform (1.4) we then also find two 
explicit formulas for the formal multiple Wilson polynomials (see Subsection 3.2). One of them 
is in terms of Kampe de Feriet series [24] . 

In Section 4 we only consider the cases where we obtain real orthogonality conditions, namely 
multiple Wilson and multiple Racah. By some limit relations we then recover hypergeometric 
representations for the examples of multiple orthogonal polynomials of type II, mentioned above. 
We also introduce some new examples like multiple dual Hahn, multiple continuous dual Hahn, 
and multiple Meixner-Pollaczek. As a result, we finally construct a (still incomplete) multiple 
AT-Askey table similar to the well known Askey scheme for classical orthogonal polynomials. 

2 Jacobi and Wilson polynomials 

2.1 Formal (shifted) Jacobi polynomials 

The (shifted) Jacobi polynomials P^ 1 '^ are a classical example of continuous orthogonal poly- 
nomials. Suppose a, (3 > —1, then these polynomials are orthogonal with respect to the Jacobi 
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weight function w <y(x ^\x) = x a (l — x)^ on the interval [0, 1]. These polynomials have the explicit 
expressions [9] 

P^\x) = ij ^*Fi(~ n,a l P + + i n + 1 U), (2-1) 
= ^±^(l-x)^ 2 F 1 f a + 1 + n '- /3 - n x], (2.2) 

where the second expression is obtained by Euler's formula [1, 15. 3. 3], [15]. We claim that the 
(shifted) Jacobi polynomials P^'^ are still formal orthogonal polynomials if we allow complex 
parameters a, (3, a + /? + lG<C\ { — 1, —2, . . .} in formula (2.1). This was already mentioned in 
[19, Theorem 2.1], but we give a different proof. 

In order to prove this claim, we require an integral representation for the meromorphic 
continuation in both variables of the beta function. Recall, e.g., from [15, § 1.1] that T is 
meromorphic in C, with simple poles at 0, —1, —2, and hence the beta function [15, § 1.5] 

. T(z)T(w) 
B(z,w) = ; v ' 
y J T(z + w) 

is meromorphic both in z and w, with simple poles at z,w = 0, —1, —2, .... From [15, § 1.5] we 
have the integral representation 



B{z,w)= f t z - l {l-t) w - l dt, »(«)> 0, »(«;)> 0. 
J o 



(2.3) 



In order to obtain a representation valid for general z, w G C\Z (compare with the Pochhammer 
formula [15, 1.6.(7)]), we follow [16, § 3.4] and consider three sheets SI, S2 and S3 of the 
appropriate Riemann surface for the function i 2_1 (l — t) w ^ 1 (in the variable t) so that 

-7r < arg(t) < 7r, -7r < arg(l - t) < n, for t & Sl \ {(— oo, 0] U [1, +oo)} , 

< arg(t) < 2vr, vr < arg(l - t) < 3tt, for t G 52 \ {[0, +oo) U [1, +oo)} , 

vr < arg(t) < 3vr, < arg(l - t) < 2tt, for t G 53 \ {(-oo, 0] U (-oo, 1]} . 

Furthermore we choose a closed contour S as in Figure 1 where a\ , (j 2 , 03 are the transition 
points between the three sheets. Note that the function t z ~ 1 (l — t) w ~ 1 is analytic on S. For the 
beta function we then have 

B{z,w) = (1 - e 2 ™)-\\ - e 2 ™)- 1 J t* -1 (l - i) 10 " 1 *, z,weC\Z. (2.4) 

Indeed, if 3ft(z) > 0, K(u;) > then the path of integration in (2.4) can be shrinked in order to 
approach the interval [0, 1], leading to formula (2.3). In particular, if z and/or w is a strictly 
positive integer, we can obtain B(z,w) by taking limits in (2.4). 

Now we prove our claim that (shifted) Jacobi polynomials with complex parameters are 
formal orthogonal polynomials. 

Theorem 2.1 Let a, (3, a + [3 + 1 G C\{ — 1, —2, . . .}, and consider the (complex) measure fi^ a '^ 
defined by 

[ h{x)d^\x) = lim (1 - e 27viz )-\l - e 2 ™)" 1 f h{t)t z {\ - t) w dt. 

Then forms a perfect system, with the corresponding nth formal orthogonal polynomial 

given by the (shifted) Jacobi polynomial P^ 1 '^ . 
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Figure 1: The contour S on the appropriate Riemann surface for the function t z (1 — t 



Proof. From (2.4) we obtain for the fcth moment 



/ 



x k dfiW>(x) = B(a + l + k,f3 + l). 



The restrictions on the parameters a and (5 guarantee that the determinant of the moment 
matrix, 

det(i?(a + r + *- 1,0 + 1)) = ff If + f °* + °, TT (--r), 



8=1 



Kr<s<n 



is different from zero. A proof for this formula uses Theorem 1 and 2 in [10, 11] (the moment 
matrix in the non-shifted case can be found in [25, 6.71.5]). Thus y^ a '^ forms a perfect system. 
Moreover, for < j < n — 1 we have according to (2.1) and (2.4) 

{l-x) j P^ l3) {x)d^\x) 

(a + l) n ^ (-w)fe(q + (3 + n + l) fc . 

2^ ~ , Tv □ B(a + A; + l,/3 + j+l) 



??■! 



fc=0 



(a + l)fcfc! 



(-l) n r(a + n + l)r(/3 + j + l) ^ (-1 



T(a + 13 + n + l) 



\TL—k 



k\(n-k)\ 



\\ {a + (3 + k + £). 



The sum on the right hand side is the divided difference Qj^n 

[0, 1, n] of the polynomial gj, n (x) = 
T\?=j+2( a + (3 + x + 1) of degree n — j — 1. Since n — j — 1 < n, this is equal to 0, which proves 
that the Jacobi polynomial is a formal orthogonal polynomial with respect to ^ a '^ . □ 



2.2 Formal Wilson polynomials 

In [27] J.A. Wilson introduced the (formal) Wilson polynomials 



p n (t 2 ;a,b,c,d) = (a+b) n (a+c) n (a+d) n 4 F 3 



-n, a + b + c + d-\-n — 1, a — t, a + t 
a + b,a + c,a + d 



1 . (2.5) 



Due to Whipple's identities [2, Theorem 3.3.3], one can show that these formal Wilson polyno- 
mials are symmetric in all four complex parameters a, b, c, d. Furthermore, with some conditions 
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on these four parameters, the polynomials satisfy a complex orthogonality with respect to the 
Wilson weight function 

Fja + t)T(a - t)T(b + t)T(b - t)T(c + t)T(c - t)T(d + t)T(d - t) 
w[t,a,o,c,a) - r(2t)T(-2t) 

Suppose that 

2a, a + b, a + c, a + d, 2b, b + c, b + d, 2c, c + d,2d(£ {0, -1, -2, . . .}, (2.6) 

so that the Wilson weight has only simple poles, and that 

a+b+c+d^ {0,-1,-2,...}. (2.7) 

Furthermore denote by the contour C the imaginary axis deformed so as to separate the increasing 
sequences of poles ({a + k}^L Q ,{b + k}^L Q ,{c + k}^L Q ,{d + k}^L Q ) from the decreasing ones 
({-a - k}f =0 , {-b - k}f =0 , {-c - k}f =a , {-d - k}%L Q ), and define the Wilson measure ^ b > c > d ) 
by 

J h{t) d^ a > b > c ' d \t) = J h(t 2 )w(t;a,b,c,d)dt. (2.8) 
Wilson [27] shows the complex orthogonality relations 

p m (t; a, b, c, d)p n {t; a, b, c, d) d^ a ^ d) (t) = 5 mn 2iM n , 



where 

M n = 2im\ (a + b + c + d + n-l) n 

r(a + b + n)r(a + c + n)T(a + d + n)T(b + c + n)T(b + d + n)T(c + d + n) 

T(a + b + c + d + 2n) 

The perfectness of the singleton system ^( a > 6 > c ' rf ) follows from Lemma 3.4 below. 

In some cases we obtain real orthogonality conditions with respect to positive measures on the 
real line [27]. If K(a), K(6), K(c), 5R(gT) > and a, b, c and d are real except possibly for conjugate 
pairs, then C can be taken to be the imaginary axis and we obtain the real orthogonality 

2 



poo 

/ p m (-x 2 ;a,b,c,d)p n (-x 2 ;a,b,c,d) 
Jo 



T(a + ix)T(b + ix)T(c + ix)T(d + ix) 



T(2ix) 



dx = 5 mn M n . 



Another case is when a < and a + b,a + c,a + d are positive or a pair of complex conjugates 
occurs with positive real parts. We then get the same positive continuous weight function 
where some positive point masses are added. In these cases we obtain the Wilson polynomials 
W n (x 2 ; a, b, c, d) := p n {—x 2 ; a, b, c, d) (see, e.g., [17]), which are real for real x 2 . 

Finally in the case that one of a + b, a + c, a + d is equal to — N + e, with N a nonnegative 
integer, we obtain a purely discrete orthogonality after dividing by T(—N + e) and letting e — > 0. 
If we take the substitution t — > x + a and the change of variables a = a + b—1,/3 = c+ d— 1, 7 = 
a + d — 1,5 = a — d we then obtain the Racah polynomials 



R n (X(x);a,f3,j,5) = 4-F3 



— n, n + a + (3 + 1, — x, x + 7 + 5 + 1 
a+ 1,0 + 6+ 1,7 + 1 



where X(x) = x(x + 7 + 5+1) and a + 1 = — N or (3 + 5 + 1 = —N or 7 + 1 = — N. They satisfy 
the discrete orthogonality 

f (a+lUl + lUP + S+lUl + S+l) x ((l + S + 3)/2) x (Mx y a g J§)(Mx)) k_ 
^ {-a + 7 + 5 + l) x (-/3 + 7 + 1M( 7 + 5 + l)/2) x (5 + l) x x\ Kn{X{Xh °' ^' 7 ' " j (A(X) j " °' 
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k = 0, . . . , n — 1. Necessary and sufficient conditions for the positivity of the weights are quite 
messy. An example of sufficient conditions is given in [27, (3.5)]. 

The formal Wilson polynomials contain as limiting cases several families of orthogonal poly- 
nomials [27, Section 4] [17, Chapter 2] like Hahn, dual Hahn, Meixner, Krawtchouk, Charlier, 
continuous Hahn, continuous dual Hahn, Meixner-Pollaczek, Jacobi, Laguerre and Hermite poly- 
nomials. Together they form the Askey scheme of hypergeometric orthogonal polynomials. 
Moreover, there exist ^-analogues [17] for all these polynomials. The results of this paper should 
probably extend to this more general case. 



2.3 Formal Wilson polynomials as a Jacobi transform 



We now recall an integral relation between the Jacobi and the formal Wilson polynomials. We 
also give a short proof which will help us to find explicit expressions for the formal multiple 
Wilson polynomials in the next section. 

Theorem 2.2 (T.H. Koornwinder) Suppose that the conditions (2.6) and (2.7) hold, and 
that 5ft(c + d) > 0, < |»(*)| < K(a). Then we have 



p n (t 2 ;a,b,c,d) = k n I P^{u)w^(u)K(u,t;a,b,c,d) du, 
Jo 



(2.9) 



with a = a + b — 1 and (3 = c + d — 1, the Jacobi weight function w^ a, ^\u) = u a (l — , the 
constant k n = n\ T(a + c + n)T(a + d + n) and the kernel (1.5). 

Remark 2.3 T.H. Koornwinder already mentioned this Jacobi transform for the Wilson poly- 
nomials in [18, (3.3)]. In order to see that these formulas coincide, one has to make some 
parameter changes and substitute 1 — th 2 s — ► u. Notice also that (2.9) is a particular case of a 
formula due to Meijer [21, p. 103]. 

Proof. By comparing the explicit formulas (2.1) and (2.5) we see that it is sufficient to prove 
that, for ££N, 



f 

Jo 



e„Ja-l,c+d-l) 



U W 



(u)K(u, t; a, 0, c, d) du 



(a - t)i(a + t) e 



(2.10) 



T(a + c + £)T{a + d + £)' 

By definition of the kernel (1.5) we have 

K(u, t; a, 0, c, d) = (a — t)i(a + t)e K(u, t;a + £, 0, c, d), 

and u i w ( - a ~ 1 ' c+d ~ 1 ^ (u) = w <ya+ ^ 1,c+d ^ 1 \u). Furthermore, Euler's formula gives the symmetry 
K{u, — t; a + £, 0, c, d) = K(u, t;a + £, 0, c, d). By replacing a + £ by a and possibly t by —t, we 
see that it remains to prove that, for < R(i) < K(a) and 3?(c + d) > 0, 



r w^+t-V ( u )K(u, t; a, 0, c, d) du = - 1 
Jo T(a + c)T(a + d) 

Denoting the left hand side of (2.11) by J, we have by definition of the kernel 



(2.11) 



J = 



w 



(a-t-l,c+d-l) 



(u) 2 Fl 



c — t,d — t 
c + d 



1 — u I du. 



T(a-t)T(a + t)T(c + d) 
In order to exchange the order of summation and integration, we notice that, if R(C— A— B) > 0, 



lim max 

s/-»i- «e[o,i] 



2^1 



A,B 
C 



yv - 2 Fi 



A,B 
C 



= 0, 
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which follows by observing that rtc+fcAi+fc) = ^ +jB_C_1 (1 + C*^ 1 )), k-oo, by Stirling's 
formula. Consequently, using the assumptions K(2i) > 0, 3? (a — t) > 0, 5R(c + eZ) > together 
with (2.3) we obtain by uniform convergence 



c — t,d — t 



T(a-t)T(a + t)T(c + d)y^i-J V^i^ c + d 

= lim V - {C ~ t)k{d J )k j M; , f w^-^ k ^(u) du 
r(a - t)r(a + t)r(c + d + k)k\ J V ; 

c — t,d — t 



y(l — u) ] du 



= lim 2 -Pi 

y-i- r(a + i)r(a + c + d - t) \a + c + d-t 

By assumption on the parameters, 9i(a + c + ti — t — (c — t + <i — t)) > 0, and hence we get from 
the Gauss formula [2, Theorem 2.2.2], [1, 15.1.20] 

c — t,d — t 



~ T(a + t)T(a + c + d-t) 2 1 \a+c+d-t 



F(a + d)T(a + c) ' 

as claimed in (2.11). □ 

3 Multiple Wilson as a Jacobi-Pineiro transform 

3.1 Jacobi-Pineiro with complex parameters 

The Jacobi-Pineiro polynomials are defined by a Rodrigues formula 

where n\ = Ilj=i n j'- ^ is wen known [23] that, provided that ocj > —1, (3 > — 1 and — Z 
for £ ^ k, the Jacobi-Pineiro polynomials are orthogonal with respect to the (positive) Jacobi 
weights ty( a i' /3 ), j = 1, . . . , m, on the interval [0, 1]. Notice that these weights form an AT system, 
and hence we obtain a perfect system of measures. Similar to Theorem 2.1, we show below that 
for complex parameters we keep formal multiple orthogonal polynomials of type II. Here we use 
the measures yS- a i'^ of Theorem 2.1 which have as support the contour S, but can be reduced 
to complex Jacobi weights w^ a ^^ on the interval [0, 1] in the case $t(otj) > — 1, §£(/?) > —1. 

Theorem 3.1 Let a>j,/3,aj + (3 G C \ {—1,-2,...}. Then the measures ^ a i^\ j = l,...,m, of 
Theorem 2.1 form a perfect system. The corresponding formal multiple orthogonal polynomial 
of type II with respect to the multi-index n is given by (3.1). 

Proof. From Theorem 1 and 2 in [10, 11] we obtain for the determinant of the matrix of 
moments (1.2) the expression 

(\n\ \ ( m nj 



X 



n n (-»■) n nn>-« 

\j=l l<r<s<rij I y l<i<j<m s=l r=l 



For our choice of parameters, this expression is different from zero, and hence every multi-index 
is normal. 
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Our claim on the (formal) orthogonality of Jacobi-Piheiro polynomials will be shown by 
induction on m. For m = 1, equation (3.1) reduces to (2.1), and the claim follows from Theo- 
rem 2.1. For m > 1, we observe first that, by the Rodrigues formula (3.1), 

P @,P), t) = t-^jl-tyf 3 cP 1 / +n _ t)fS+nip (( a2 ,..., a T ),/3+m) ( ^ , 3 2) 

n w dt n i V v ; (ri 2 ,...,n m ) v > 

By induction hypothesis, p^ 2 >--^ Q ™)'^+ ni ) j s a polynomial of degree \n\ — m, and thus there 
exist scalars Cj with 

|n|-ni 

p((a 2 ,...,Q!m),/3+ni)/,x _ p(ai+ni,/3+ni)/,A 
Un 2 ,...,n m ) W ~ C i^j W- 

j=0 

From the Rodrigues formula for to = 1 and (3.2) we conclude that 

|n|-ni 

i=o 

implying that pj&'W is a polynomial of degree |n| with J t k P^'^ {t) d^ ai) ^ = for = 
0,1, ...,ni — 1 by Theorem 2.1. The other orthogonality conditions are obtained by observing 
that (3.1) remains invariant if one changes the order in the product in (3.1). □ 

With help of the Leibniz rule applied to the Rodrigues formula, the authors in [14] derive 
for to = 2 an explicit expression in terms of m = 2 sums. We now give a generalization of this 
formula for m > 2, using the notation 



M p;r : 

q,n 



/; 91 ■ ■ ■ ■ ■ 9r 

<i>\i>\ •■ ■ ■ ■ •■ 



V 



x (3.3) 



"l n m II (ft)\%\ I! (5 , «,l)|fc|_fc 1 ' ' ' (9i,in~l)k m m kj 

'■= Yl ••• E ¥ — U(~ n jhfj, 

k 1=0 ^ fc m =Q FJ fj (A,l)\k\- kl " " " (^i,m-l)k m i= l 3 ' 

to sums 1-1 % ~ 

where k = (ki,...,k m ), n G (NU{0}) m , / € C p , ^eC« and G C" 1 " 1 . 

We also give in (3.5) another new explicit expression for the formal Jacobi-Piheiro polynomials 
which reduces to (2.2) if m = 1. 

Theorem 3.2 Zei e = (1, . . . , 1) &e a multi-index of length m and s(n) = (ni, n\ + rc-2, . . . , |n|). 
Denote by v^^ the vector v without the jth component. For the J acobi- Pineiro polynomials we 
have the hypergeometric representations 



p (S,/3) m = (a + e)n M i;2 ( («i + + ni + 1); (a + n + e)( m ) : (a + s(n) + ((3 + 
™ l J n! ( ai + i) ; (5 + e)(i) : (a + s (n) + (/3 + l)e)(™) 



tej 
(3.4) 



anc? 

P^(t) = (1 " ^ m+ ,F m ( d + H + c ; " I*' t) , (3.5) 

w/iere n! = flyli n j ] - an d (« + e)n = n7Li( a j + l)rij • 
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Proof. We prove (3.4) and (3.5) by induction on m. For m = 1, equations (3.4) and (3.5) 
reduce to (2.1) and (2.2), respectively. 

In case m > 2, we use formula (3.2), where the induction hypothesis enables us to express 

the right hand polynomial P^ ^ +ni ^ as a hypergeometric sum. After exchanging the order of 
summation and differentiation (for \t\ < 1 in case of formula (3.5)), we apply the formulas 

t"«i(l _ t)~ ^ ^i+-i+|fe (1) l(l _ t f+n^ (3. 6 ) 

^ ("i + g + "i + ("i + "i + i)|g|- fcl (- ni ) kl S 

and 



t - ai d^ t a 1+ n l+ i = (gl + l)»i(«l +"1 + l)i t i 

df 11 (ai + l)» 



(3.7) 



and obtain the right-hand expressions of (3.4), and (3.5), respectively. It remains to prove the 
claims (3.6) and (3.7), the second one being obvious. We observe that the left-hand side of (3.6) 
can be transformed using the Rodrigues formula for m = 1 and (2.1), leading to the expression 

= K + iff') | + iU t ( " ni)fcl (ai + ' + P + ni + t\ 

^ («i + /3 + nx + l) |g| ( ai + ni + l) |g| _ fci (-m) fel #l 

as claimed in (3.6). □ 

In the above proof we have shown implicitly that the right hand side of (3.5) is a polynomial of 
degree < \n\ in t. 

With the explicit formula (3.4) it is possible to compute the recurrence coefficients by com- 
paring the highest coefficients in the recurrence relation of order m + 1, see [14] for m = 2. The 
resulting formulas are quite complicated, and thus we omit further details. 

3.2 Formal multiple Wilson polynomials 

We now consider m Wilson weights 

w(-;a,bj,c,d), j = l,...,m, b E -b k gZ, £ / k, (3.8) 

defined as in (1.3), that is, we change only one parameter (recall the symmetry of the Wilson 
weights in all four parameters). Similar as in the scalar case (Wilson polynomials) we suppose 
that for j = 1, . . . , m 

2a,a + b j ,a + c,a + d,2b j ,b j + c,bj + d, 2c, c + d, 2d {0,-1,-2,...}, (3.9) 

so that the m Wilson weights have only simple poles, and that 

a + bj + c + d^ {0,-1,-2, ...}. (3.10) 

As in the scalar case, we write yS- a,bj,c,d ^ for the resulting measures, where it is possible to choose 
a joint contour C which is the imaginary axis deformed so as to separate the increasing sequences 
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of poles ({a + k}%L , {b x + k}™ =0 , ...,{b m + k}™ =0 , {c + k}™ =0 , {d + k}%L ) from the decreasing 
ones ({-a - fc}~ , {-fc - , . . . , {-6 m - fc}~ , {-c - {-d - fc}^o). 

Let us show that the (possibly complex) Wilson measures ^,( a ' b i' c ' d ) form a perfect system. 
The corresponding multiple orthogonal polynomials will then be referred to as formal multiple 
Wilson polynomials. A basic observation in what follows is that, under some additional con- 
ditions, the formal multiple Wilson polynomials can be written as a Jacobi-Piheiro transform, 
similar to (2.9). 

Theorem 3.3 Suppose that (3.9) and (3.10) hold and that bg — b^^Llj,l^k. Then the Wilson 
measures ^,( a ' b i' c ' d ) for j = 1, . . . , m form a perfect system. 

7/K(a) > and 3?(c + d) > 0, the formal multiple Wilson polynomial with multi-index n can 
be written as 

I' 1 

p H (t 2 ;a,b,c,d) =k H pf' P] ^(u)w {a ~ 1,/3) (u)K(u,t;a,0,c,d) du, (3.11) 
Jo 

for < |K(t)| < 5ft(a), where a = (a + b\ — 1, . . . , a + b m — 1) = (a — l)e + b and (3 = c + d — 1. 
The normalizing constant k^ = n\ T(a + c + |n|)r(a + d + |n|) is chosen so that it corresponds 
with (2.9) in the case m = 1 and the kernel K(u,t;a,b,c,d) is defined as in (1.5). 

Before we prove this theorem we need some technical lemmas. 

Lemma 3.4 The system ofm measures fii, . . . , fi m is perfect if and only if, for every multi-index 
n, there exists a polynomial p^ of exactly degree \n\ so that 

J ' Pn(z)z k dfij(z) = 0, k = 0,...,nj-l, j = l,...,m, (3.12) 

and 

J PH (z)z n i d H (z) + 0, j = 1, . . . , m. (3.13) 

In this case, pa is the (up to normalization unique) multiple orthogonal polynomial of type II 
with respect to n. 

Proof. Suppose first that fi m is perfect, and take as p^ the multiple orthogonal polyno- 

mial of type II with respect to n. Then it only remains to verify (3.13). Suppose the contrary, 
that is, J pfi(z)z nj dfij(z) = for some j. Then p^ is also a multiple orthogonal polynomial of 
type II with respect to V = (ni, rij-i, rij + 1, nj+i, ...,n m ), in contrast to the normality of the 
multi-index V. 

We will prove the other implication of this lemma by showing by induction on the length 
|n| that n is normal. The multi- index is always normal, suppose therefore that n is of length 
> 1, with its jth component strictly greater than 0. Let be a multiple orthogonal polynomial 
for n. If q^ would have degree strictly less than \n\, then it would be a multiple orthogonal 
polynomial for the multi-index V = (m, ...,nj-i,rij — 1, n^+i, n m ) with length |n| — 1. By 
normality of V we then have that there exists a nonzero constant c so that q^ = cpp. Thus 
J qfi{z)z n i~ l dfXj(z) by (3.13), in contradiction to the orthogonality relation (3.12) for q^. 
As a consequence, q^ has the precise degree \n\, and thus ft is normal. □ 

Lemma 3.5 Suppose that the singleton systems form a perfect system for j = 1, ...,m, with 
corresponding orthogonal polynomials {pn^} n - Then the system of m measures fi\, . . . , fi m is 
perfect if and only if, for every multi-index n, there exists a polynomial Pfi and scalars c^ k so 
that 

\n\ 

Mx) = E <%Pi\ x )> C S^0, cgj^O, j = l,...,m. (3.14) 
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In this case, pa is the (up to normalization unique) multiple orthogonal polynomial of type II 
with respect to n. 

Proof. If n is normal and Pft is the corresponding multiple orthogonal polynomial, then (3.14) 
follows by taking 



where we observe that the denominator is nontrivial according to (3. 12), (3. 13) for the singleton 
system fij. In addition, c^j, = for t < rij by (3.12), c^. / by (3.13), and cJL / 0. 

Conversely, (3.14) plus the perfectness of the singleton system fij implies (3. 12), (3. 13), and 
hence the system fii, . . . , fj, m is perfect by Lemma 3.4. □ 

We now prove Theorem 3.3 by showing that (the analytic extension of) the integral expression 
(3.11) is a possible candidate for a formal multiple Wilson polynomial. 

Proof of Theorem 3.3. According to the assumptions (3.9) and (3.10) of Theorem 3.3, we 
find that ctj, (3, ctj + (3 + 1 G C \ { — 1, —2, . . .}, j = 1, . . . , m, and at — a& Z whenever £ ^ k. It 
follows from Theorem 3.1 that the Jacobi-Piheiro system ^- a ^^\ j = 1, ...,m, is perfect. From 
Lemma 3.5 we may conclude that there exists scalars ci k so that 

pf^\t) = £ cgp^V), 4£ 3 / 0, eg},, ^ 0, j = 1, . . . ,m. 

From, e.g., (3.4) we see that Jacobi-Piheiro polynomials are rational functions in each of the 
parameters ctj or (3. Taking into account (2.4) and (3.15), we may conclude that any of the 

coefficients is a meromorphic function in each of the parameters ctj or (3. We now introduce 

q H (t 2 ;a,b,c,d) = k H pf ,/3) '(u)^ " 1 ^ ] (u)K(u,t;a,0,c,d) du. 

Jo 

This function is well defined for < |5ft(i)| < 5R(a) if 5R(a) > and 3ft(c + d) > 0. However, using 
(2.9), we obtain for every j = 1, . . . , m that 

1*1 ,i 

fe (t 2 ; a, b, c, d) = kft V 4 J J / P^' /3) (x)u;^'' 3 )(n)^(n, t; a, 6„ c, a!) du, 

l=n 3 ' ^ 

Inl 

- V ni na + c+\n\)T(a + d+\n\) U) 2 



and thus 



\n\ 



qn(t 2 ;a,b,c,d) = ^2 -ft( a + c + t)\fi\-e (a + d + l)\fi\-i cl J ] p e (t 2 ;a,bj,c,d) 
e= nj 

\n\ 

= ^2d { ^p e (t 2 ;a,b 3 ,c,d). (3.16) 



Observing that the expressions on the right-hand side of (3.16) are polynomials in t and mero- 
morphic in any of the parameters a, 6, c, d, we see that the right hand expression of (3.16) is well 
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defined and independent of j for any choice of t and of the parameters a, b, c, d, as long as (3.9) 

l n,\n\ 



and (3.10) hold and b( — Z, £ / A;. Moreover, the new coefficients and for d^L, are 



different from zero. 

Thus, qfi(t 2 ;a,b,c,d) defined by (3.16) is a suitable candidate for a formal multiple Wilson 
polynomial, and the system of Wilson measures is perfect by Lemma 3.5. □ 

We now want to deduce explicit expressions for the formal multiple Wilson polynomials 
based on the explicit expressions (3.4) and (3.5) for the Jacobi-Piheiro polynomials. Here we 
use the Kampe de Feriet series [24] 



tpp-pi;P2 
9:91:92 



/ : 9; h 



V PI P2 

00 n (/<)* k n (9e)k~j n (^)i 



x,y 



E e=i \ ^ e=i 
9 2-^i 91 



1=1 



x k 3 y 3 



n (<Pdk 3=0 n (^)fc-j ft mi {k 3)1 Jl ' 

i=\ 1=1 e=i 



(3.17) 



which are a generalization of the 4 Appell series in 2 variables. Notice that, for p = q = 0, the 
Kampe de Feriet series is a product of two hypergeometric series. Also, in the case m = 2, our 
functions MF^- defined in (3.3) are (finite) Kampe de Feriet series 



M p ' r ; , 



f; gi 9r 
4>\ fa ■ ■ ■ ■ ■ A 



(xi,x 2 ) 



rpp:l;r+l 
r q:0;r 



f ■ (~ni); (-n 2 ,gi,...,g r 

$ ■ 0; (-01, -,Yv) 



X]_,X 2 



(3.18) 



In what follows the parameters in the Kampe de Feriet series will always be chosen so that the 
sum in (3.17) is finite, and hence we are not concerned with convergence problems. 

Corollary 3.6 Let e = (1, . . . , 1) be a multi-index of length m, s(n) = (ni,ni + n 2 , ... ,\n\) and 



a + bj + c + d — 1, j = 1, . . . , m. Denote by the vector v without the jth component. 

(3.19) 



For the multiple Wilson polynomials we have the hypergeometric representations 
Pn(t 2 ; a, b, c, d) = (ae + b) n (a + c)\ H \ (a + d)\ H \ 



x M 



3;2 
3,n 



(a-t,a + t,a! + n x ); (ae + b + n) (m) : (a + s(n))^ 
(a + c,a + d,a + 6i);(ae + 6)( 1 ) : (<r + s(n))( m ) 



and 



Pn(^ 2 ; a, &, c, d) = (ae + 6) fl (a + c)\ H \ (a + d)|^| 



^2:i ; m+i | (a — i, a + i) : (c + d — 1) ; (ae + 6 + n, 1 - c - d - \ft\ 
2:0;m I (a + c,a + d) : ();(oe + S) 



(3.20) 
!.d 



w/iere (ae + 6)^ = n^i( a + ^0% 



Proof. First note that (1 — x) & = Yle=o(P)ljri wn i cri converges in the unit disk, so that the 
expression (3.5) can then be written as 



P. 



{5,13) 



(a + e) H :i;m+i /" ~ : (/?); (« + ™ + e, -P - \n\) 

=T -Fn 



n! 



0:0;m 



: -; (a + e) 



(3.21) 



Then start from the Jacobi-Piheiro transform (3.11) and replace the Jacobi-Piheiro polynomial 
pC"-/ 3 ) ky eX pii c it expressions (3.4) and (3.21), respectively. Since the sums are finite, we can 
interchange the integral with the sums. Applying (2.10) then completes the proof. □ 



For m = 2, we may apply (3.18) to (3.19), leading to a representation as a Kampe de Feriet 

vial 1 

(3.20) as a Kampe de Feriet series of type i^-o-'l ■ 



3" 1 '3 

series of type F 3 .' .' 2 ■ It seems to be non-trivial to derive from this formula the representation 
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4 Limit relations 



In this section we consider some cases in which the orthogonality conditions of the formal multiple 
Wilson polynomials reduce to orthogonality conditions with respect to a positive measure on the 
real line. We then recover the multiple Wilson and multiple Racah polynomials. Next we use the 
limit relations between the orthogonal polynomials in the Askey table [17] to obtain some new 
examples of multiple orthogonal polynomials and some known examples. In particular we look 
at what happens with the explicit expressions (3.19) and (3.20) applying these limit relations. 
Most of these examples are known to be AT systems (which implies that every multi-index is 
normal) . 



4.1 Multiple Wilson 

With some restrictions on the parameters the orthogonality conditions of the formal multiple 
Wilson polynomials reduce to real orthogonality conditions with respect to positive measures 
on the real line. Let bj > 0, j = 1, . . . , m, b% — ^ Z whenever I ^ k, 3?(a), 3i(c), 3?(d) > and 
a, c, d be real except for a conjugate pair. In that case the imaginary axis can be taken as the 
contour C. The multiple Wilson polynomials 

W H (x 2 ; a, b, c, d) := pn{-x 2 ; a, b, c, d) (4.1) 

then satisfy the real orthogonality relations 

/>oo 

/ (x 2 ) k W a (x 2 ;a,b,c,d) 
Jo 

k = 0, . . . , rij — 1, j = 1, . . . , m. If a < 0, a + bj > 0, j = 1, . . . , m, and a + c, a + d are positive 
or a pair of complex conjugates with positive real parts, then we obtain the same orthogonality 
conditions but with some extra positive point masses. 



r(a + ix)T{bj + ix)T(c + ix)T(d + ix) 
T(2ix) 



dx = 0, 



4.2 Multiple Racah 

As in the scalar case it is also possible to obtain a purely discrete orthogonality. The multiple 
Racah polynomials Ra{-; a, (3, 7, 5) (where we only change the parameter a with ct£ — at Z 
whenever t 7^ k) satisfy the discrete orthogonality 

f K + 1)*(7 + !)«(/? + * + 1)«(7 + S + M7 + 8 + 3)/2) a R , x(x) . a §)(x(x))k _ Q 
to + 7 + ^ + l),(-/3 + 7 + 1),((7 + 5 + l)/2) x (5 + l) x x\ n ^ ^ ' ^ 7. Warn - U, 

A; = 0, . . . , n — 1, j = 1, . . . , m, where 

A(x) =x(x + 7 + (5 + l) and fi + 8+1 = -N or 7 + 1 = -N. 

They can be found from the polynomials Pft(t 2 ; a, b, c, d)/((ae + 6)^(a + c)|^| (a + by the 

substitution t —>■ x + a and the change of variables aj = a + bj — 1,(3 = c + d — 1,7 = 
a + ci — l,<5 = a — ci. For the multiple Racah polynomials we then have the expressions 



3,n 



_ r 2:0;m 



Rft(\(x);a,(3,j,5) 

(-x, x + j + 6+l,a 1 + P + n 1 + l);(a + n + e) {m) : (a + s(n) + {(3 + l)e)W 
(/3 + <5 + 1, 7 + 1, ai + 1); (a + e)^ : (a + s(n) + (/? + l)e)( m ) 

,n2:i;m+i ^ (~x, x + 7 + 5 + 1) : (a + n + e, -/3 - |n| 
(/? + ,5+l,7 + l):0;(a + e) 

In this case, necessary and sufficient conditions on the parameters to have positive weights are 
quite messy. 



1,1 
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Remark 4.1 Remember that the Wilson weight is symmetric in the four parameters so that 
we can switch these parameters in the change of variables. We then obtain multiple Racah 
polynomials (of type II) where we change other parameters. However, all these cases are related 
to each other. For example we have multiple Racah polynomials Rfi(-; a, f3, 7, 5) where we only 
change the parameter (3 in the weights (Pi — Pk Z whenever i 7^ k). In that case we have that 
a + 1 = — N or 7 + 1 = — N. As a second example it is possible to change the parameters /3, 7 
and 5 in such a way that Sj + jj and 5j + (3j don't change and 7^ — 7fc Z whenever t 7^ k. 
Here we assume that a + 1 = — N or f3j + 5j + 1 = — N. We then denote these multiple Racah 
polynomials by R n (-; a, (3, 7, 5). However, we don't find another family of polynomials because 



and 



Rfi(X(x); a, (3, 7, 5) = R H (X(x); (3 + 5e,a- 5, 7, 5) 



R H (X(x);a,P,j,S) = R H (X(x);^,a + f3j -7^,0, 7^ +Sj - a). 



(4.2) 



(4.3) 



These relations will help us in some of the examples of the subsections below to find explicit 
expressions for the polynomials. 



4.3 Some new examples 

4.3.1 Multiple continuous dual Hahn 

Let hi — bk Z whenever I 7^ k. The multiple continuous dual Hahn polynomials satisfy 
the orthogonality conditions of the multiple Wilson polynomials where we let d — ► +00 (after 
dividing by T(d) 2 ). Similarly we obtain real orthogonality conditions with respect to a positive 
measure if bj > and a, c are positive or a pair of complex conjugates with positive real parts. 
We then denote the nth multiple continuous dual Hahn polynomial by Sfi(-;a,b,c). These 
polynomials satisfy the orthogonality conditions 



f 

Jo 



(x 2 ) k S H (x z ;a,b,c) 



T(a + ix)Y(bj + ix)T(c + ix) 



k = 0, . . . , rij 

Sfi(x 2 ; a, b, c) = lim 



T(2ix) 

1, j = 1, . . . , m. It is clear that 
Wft(x 2 ;a,b,c, d) 



dx = 0, 



(a + rf)|fj| 

so that the multiple continuous dual Hahn polynomials have the explicit expressions 

(a — ix, a + ix); (ae + b + n)( m ) 



(4.4) 



Pn{x ;a,b,c) 



{ae + b)fi{a + c)\ H \ M%1 



(a + c,a + 6i);(ag , + S)( 1 ) 



( -,iS l 1 \ Z72:0;m I (d ~ IX , (1 + ix) I (); (ae + b + n) 

(ae + b) H {a + c) N F 1:0m ^ 

V (a + c) : ();(ae + 6) 



1,-1 



4.3.2 Multiple dual Hahn 

Consider 7j,<5j, j = 1, . . . ,m, so that 7j,<5j > —1 or 7j,<5j < — TV for each j and that jj + <5j 
is independent of j. Suppose also that — "fk Z whenever I ^ k. The multiple dual Hahn 
polynomials, denoted by Rfi(-;j,S,N), satisfy the system of discrete orthogonality conditions 



N 
z=0 



(2a + 7j + ^ + l)(7j + jO £ HVVV! 
■l) a: (a; + 7 i ,-+<5j + l)j V+ i(<5 J - + l) a .s! 



Rft{\{x)rf,lN) (X(x)) k = 0, 
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k = 0, . . . , n— 1, j = 1, . . . ,m, where X(x) = x(x+j+5+l). The multiple dual Hahn polynomials 
are related to the multiple Racah polynomials like 



Rn(X(xy,i,d,N) 



lim Rfi(X(x);a,-S - (N + l)e,j,5) 

lim R^(X(x);j, a — 7,- — Sj — N — 1, a, 7,- + 5j — a), 



where we use (4.3). They then have the explicit expressions 

RAMx)^6N) = ^iJ^x + ^ + S. + iy^j + n + e)^ 
n[ { h 7 ' ' ' J 2 > H { (-N, 71 + 1); (7 + 

p 2:0;m f X + Jj + <5j + 1) : (); (7 + ft + e) 



^l:0;t 



-Ar):();(7 + e) 



1,-1 



(4.5) 



4.3.3 Multiple Meixner-Pollaczek 

The multiple Meixner-Pollaczek polynomials pi X \-;(j)) are multiple orthogonal polynomials (of 
type II) associated with the system of weights e^^~^ x \T(X + ix)\ 2 , X > 0, < <pj < it, 
j = 1, . . . , m, where the (f>i, . . . , (p m are different. These weights form an AT system on the 
positive real axis [22, p. 141]. The multiple Meixner-Pollaczek polynomials satisfy the conditions 



/ x k P^>(x; ^e^-^lTiX + ix)\ 2 dx = 0, 
Jo 



k = 0, . . . , rij — 1, j = 1, . . . , m. Similar as [17, (2.3.1)] it is easy to check that 



lim 



S H ((x - t) 2 ; X + it, t cot <j>, X - it) 
(t esc (p)fi nl 



(4.6) 



where t cot 4> = (t cot (j>i, ■ ■ ■ , t cot (f> m ) and t esc 4> = it esc (/>i,...,t esc (f>m)' The multiple Meixner- 
Pollaczek polynomials then have the explicit expression 



(2%inr=i^ 



MY. 



1:0 ( (X + ix);- 



(2A); 



e — e 



where e 2l( ^ = (e , . . . , e 2i<?im ) . Here we do not have a Kampe de Feriet representation such 
as in (3.20). 



4.3.4 Formal multiple continuous Hahn 

Similar as in [17, (2.1.2)] we can use the limit relation 

D r a\ v Pn{(t + y) 2 ;a-y,b + ye,c-y,d + y) 

Pfi{t;a,b,c,a) = hm 

y—KX 



(a + c-2y)\ n \ n\ 

in order to find the formal continuous Hahn polynomials. They have the explicit expressions 

P n (t;a,b,c,d) = (ae + b) a (a + d) lal M 2 ' n \ - m f ^u m ) 

' V (a + d, a + Oi); (ae + o^ 1 - 1 : (<7 + s(n)) ( . m - ) 

= (ae + 6)^(a + c%| 



(4.7) 



x F n 



i:i;m+i ( (a + t) : (c + d - 1); (ae + 6 + n, 1 — c — d — \n\ 



l:0;m 



(o + d) : ();(ae + 6) 



1,1 
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where aj = a + bj + c + d — 1, j = 1, . . . , to. If the parameters satisfy (3.9) and (3.10) and 
bt — bk $l7L whenever i ^ k, then these polynomials satisfy the orthogonality conditions 

J Pn{t\ a, 6, c, d)F(a + t)T(bj - t)F(c + t)F(d - t)t k dt = 0, (4.8) 

k = 0, 1, . . . , rij — 1, j = 1, . . . , m, where C is a contour which is the imaginary axis deformed so as 
to separate the increasing sequences of poles ({b± + k}^L Q , . . . , {b m + fcjfclo' + ^}fc^=o) f rom t ne 
decreasing ones ({— a— k}^L Q , {— c— k}^L Q ). In the scalar case (to = 1) it is possible to obtain real 
orthogonality relations with respect to a positive measure if we suppose K(a), 9ft(c), 9ft(eZ) > 
and a = b, c = d. This is not possible in the multiple case. Therefore it would be nice if 
we could find another family of multiple Wilson polynomials in the sense that we change two 
parameters forming a pair of complex conjugates. 



4.4 Some classical discrete multiple orthogonal polynomials 

In this section we obtain hyper geometric formulas for the classical discrete examples of multiple 
orthogonal polynomials of type II, introduced in [5], which are all examples of AT systems. In 
particular we use the limit relations between these polynomials and the Racah polynomials [17]. 



4P\r 
q,n 

terms of a Kampe de Feriet series is new (if it exists) . 



Their M^- representation is already known in the cases to = 1,2. The explicit expression in 



• Multiple Hahn: These multiple orthogonal polynomials (of type II) satisfy orthogonality 
conditions with respect to to hypergeometric distributions 



N 



^ = ^—\ {N^W 5 ^ ^>-l,P>-l, 

ol\ — oik 4- {0, 1, . . . , N — 1}, £ ^ k, on the integers 0, . . . , N. They can be found from the 
multiple Racah polynomials taking 7 + 1 = — N and S — > +00, so that 

n a-,f3;N ( , ^2,2 ( (-x, a x + + m + 1) ; (a + n + e)( ro > : (a + s(n) + (/? + l)e)( 1 ) 

W#t W yw 2,nl (-iV^i + ^^a + e)^) : (a + s(n) + (/3 + l)e)M 



pi:i;m+i / :(/?); (<3 + n + e, -/? - |n| 
1:0;m V (-iY):Q;(a + e) 



1,1 



Changing only the parameter does not give another family of polynomials because of 

Q'n N ( x ) = C Q^ a,N (N — x) with C some constant (depending on n,a and (3). However, 
we will need an explicit formula in powers of x for these polynomials to reach multiple 
Meixner I and multiple Laguerre II. Using (4.2) and the limits we mentioned above, we 
find that 



o a;0;N ( , _ ^2;2 ( « + /?! + "1 + 1); + s(H) + (a + l)e)( 1 

™ W 2 '"l (-A^a + l);(/3 + s(n) + (a + l)e)( m ) 



Applying these limits to the Kampe de Feriet representation does not work. 
• Multiple Meixner I: In this case we consider m negative binomial distributions 

= 2 k\ " L <5fc ' < c i < ^ /? > °> 

fc=0 
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with all the Cj, j = 1 . . . , m, different. We get these polynomials from the multiple Hahn 



polynomials replacing a = j3 — 1, (3j 

obtain 



N 



l — Cj 



and letting N — > +00. We then 



1:0 



(-*);() 

mo 



c — e 



where 

c 



— ( ci^ 1 c m -l \ 
\ ci ' • • • ' c m y ■ 



• Multiple Meixner II: In the case of multiple Meixner II polynomials we only change 
the parameter (3 in the negative binomial distributions, so that 



k=0 



k\ 



S k , < c < 1, /3j > 0, 



with fy — f3 k Z whenever £ ^ k. Taking aj = (3j - 1, (3 = and letting N — > +cx) in 

the explicit formulas for the multiple Hahn polynomials Q% N ■> we obtain 



1;1 / (-x);(/3 + n)M 



c- 1 



(A);G9)« 

l:0:m / : 0? (0 + «) 



0:0;m 



() = (); 09) 



c- 1 1 - c 



Multiple Kravchuk: These polynomials satisfy the orthogonality conditions (1.1) with 
the m binomial distributions 



fc=0 v 7 



N-k x 
Ok 



<Pj < 1, 



where all the pj, j = 1 . . . ,m, are different. These polynomials are related to the multiple 
Hahn polynomials Q^ ,N replacing j3 = t, aj — ► j^— t and letting £ — > +oo. We then get 



* 1 j ^ 1 (-7V); () 



where I = (±, . . . , ±). 
• Multiple Charlier: In the case of multiple Charlier we consider m Poisson distributions 

00 a k 



k=0 



with all the aj, j = 1 . . . , m, different. The corresponding multiple orthogonal polynomials 
(of type II) can be found from the multiple Meixner I polynomials taking Cj = Q a | /3 and 
letting (3 — > +oo. The multiple Charlier polynomials then have the explicit expression 



C a (x] = M i;0 ( (-*); 



where i = (±,...,jJ). 
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4.5 Some classical continuous multiple orthogonal polynomials 

In this section we recall some classical continuous examples of multiple orthogonal polynomials of 
type II where the measures (or weight functions) form an AT system and obtain hyper geometric 
formulas for these polynomials. Their M.^ representation is already known in the cases m = 
1,2. The explicit expression in terms of a Kampe de Feriet series is new (if it exists). For an 
overview of these polynomials and their properties we recommend [14, 4]. 

• Jacobi-Pineiro: In Subsection 3.1 we recalled the Jacobi-Piheiro polynomials P?" 5 , 
which, in the case aj,f3 > —1, are the multiple orthogonal polynomials (of type II) with 
respect to the Jacobi weights w aj, @(x) = x Qj (l — x)@, j = 1, . . . , m, on the interval [0, 1]. 
Here ag — a& ^ Z whenever £ ^ k. Similar as in the multiple Hahn case we have that 

Pfi*^\x) = (— \)\ n \p^ ,a \\ — x). So, changing only the parameter (3 does not give another 
family of polynomials. For these polynomials we have 



(<* + %; M i;i I (a + p 1 + n 1 + 1); + s(n) + (a + l)e)« 
n\ W I (a + 1); + s{n) + (a + l)e)M 



xe 



Multiple Laguerre I: The multiple Laguerre I polynomials L°i are orthogonal on [0, oo) 
with respect to the m weights Wj{x) = x c *e~ x , where aj > —1, j = 1, . . . ,m, and at — 

Qfc Z whenever £ ^ k. They can be found from the Jacobi-Piheiro polynomials P~'^ 
substituting x — > -| and letting (3 — > oo. We then obtain the hypergeometric expressions 



T a M (*±ek M o-,i ( ();(a + n + e)( m ) 

(a + e)fj x ^a+n+e 

6 m-fm 



xe 



n! \ a + e 

Multiple Laguerre II: In this case the polynomials L^'^ have the orthogonality condi- 
tions (1.1) with respect to the weight functions Wj(x) = x a e~ C: > x , where a > —1, Cj > 0, 
j = 1, . . . , m, and all the Cj different. They can be obtained from the Jacobi-Piheiro poly- 
nomials p( a '^ by the substitutions x — > |, taking /3j = c^t and letting i — > oo. We then 
get 



r{^ ( , _ (<* + %| ■ ,0;0 ( ();() 

^ [X) ~ n! 7Ml '"V (« + !);() 



Multiple Hermite: In the multiple Hermite case we consider the multiple orthogonal 
polynomials (of type II) PL'i with respect to the weights Wj(x) = e~ x +c i x 7 j = 1, . . . , m, on 
(— oo, +oo). Here the Cj are different real numbers. These polynomials can be found from 
the Jacobi-Piheiro polynomials P^ taking aj = (3 + Cj^/j3, the substitution x — ► 
and letting (3 — ► +oo after multiplying with some constant depending on n and (3. It seems 
that we can't hnd explicit expressions for the multiple Hermite polynomials in terms of a 
MFf- series or a Kampe de Feriet series using such arguments. 

q,n 
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Figure 2: The (incomplete) multiple AT-Askey scheme 



5 Conclusion 

In Figure 2 all these and some extra limit relations are combined. This scheme is a generalization 
of the Askey scheme in the scalar case. In each of the examples of this scheme the measures 
have the same support. Although for multiple Wilson, multiple Racah, multiple continuous dual 
Hahn and multiple dual Hahn it is still an open question, we believe that they are all examples 
of AT systems which is the reason we call it the multiple AT-Askey scheme. 

This scheme doesn't contain all the possible examples of multiple orthogonal polynomials 
which reduce to the orthogonal polynomials of the Askey scheme. In [14] the authors also 
mentioned some examples of Angelesco systems (with their hyper geometric expression). It 
is also possible to change more than one parameter in the Wilson weight (maybe with some 
correlation) in order to find other examples of multiple Wilson polynomials. Then it is for 
example possible to obtain multiple continuous Hahn polynomials (with positive measures on 
the real line) using some limit relations. 
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